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Abstract 

We consider the Hamiltonian dynamics and thermodynamics of spherically 
symmetric spacetimes within a one-parameter family of five-dimensional Love- 
lock theories. We adopt boundary conditions that make every classical solu- 
tion part of a black hole exterior, with the spacelike hypersurfaces extending 
from the horizon bifurcation three-sphere to a timelike boundary with fixed 
intrinsic metric. The constraints are simplified by a Kuchaf-type canonical 
transformation, and the theory is reduced to its true dynamical degrees of free- 
dom. After quantization, the trace of the analytically continued Lorentzian 
time evolution operator is interpreted as the partition function of a thermody- 
namical canonical ensemble. Whenever the partition function is dominated by 
a Euclidean black hole solution, the entropy is given by the Lovelock analogue 
of the Bekenstein-Hawking entropy; in particular, in the low temperature limit 
the system exhibits a dominant classical solution that has no counterpart in 
Einstein's theory. The asymptotically flat space limit of the partition function 
does not exist. The results indicate qualitative robustness of the thermody- 
namics of five-dimensional Einstein theory upon the addition of a nontrivial 
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I. INTRODUCTION 



A gravitational theory whose Lagrangian density consists of multiples of lower- 
dimensional Euler densities has the property that the field equations are second order in 
the metric [|l],@. These theories, known as Lovelock theories, include Einstein's theory with 
a cosmological constant in all dimensions greater than two, and in five or more dimensions 
they provide genuine curvature squared generalizations of Einstein's theory. Among all cur- 
vature squared generalizations of Einstein's theory. Lovelock theories therefore have a special 
status in that they preserve the number of degrees of freedom: a generic curvature squared 
action yields field equations that are fourth order in the metric, containing thus more degrees 
of freedom than Einstein's theory. This has generated wide interest in Lovelock theories, 
especially in the contexts of cosmology and black hole physics p-p^|. 



The purpose of the present paper is to analyze the classical and quantum dynamics 
of spherically symmetric Lovelock gravity by the Hamiltonian methods recently developed 
by Kuchaf ||2^. These methods have previously been applied to spherically symmetric 
Einstein(-Maxwell) gravity in four dimensions pO|-p^, vacuum dilatonic gravity in two di- 
mensions p3|j25|j26[] , and to related systems [^,^]; for related discussion, see Refs. p9|-|35 



At the classical level, we wish to find a canonical transformation that introduces the mass 
parameter of the spacetime as a new canonical variable, use this transformation to simplify 
the constraints, and reduce the theory to its true dynamical degrees of freedom. At the quan- 
tum level, we wish to derive from the quantum theory a partition function that describes 
the equilibrium thermodynamics of a Lovelock black hole in the canonical ensemble. 

The issues of prime interest are twofold. First, although Lovelock theories have the same 
set of canonical variables as Einstein's theory, the Lovelock Hamiltonian is, in general, a mul- 
tivalued function of the canonical variables One anticipates that this multivaluedness 
may introduce additional features in the canonical formulation and Hamiltonian reduction, 
even though the Lovelock analogue of Birkhoff 's theorem [§] strongly suggests that the local 
considerations should differ little from those in Einstein's theory. Second, certain Lovelock 
black holes have thermodynamical properties that differ qualitatively from those of Einstein 
black holes; in particular, a Lovelock black hole can be stable against Hawking evaporation 
in asymptotically flat space This leads one to ask whether Lovelock theories might 



admit quantum thermodynamical ensembles with boundary conditions that do not give rise 
to well-defined ensembles in Einstein's theory. 

The number of possible Lovelock terms in the action increases with the dimension of the 
spacetime, and different choices for the coefficients yield qualitatively different theories. In 
this paper we shall aim for concreteness at the expense of generality. We concentrate on a 
specific one-parameter family of Lovelock theories exhibiting both a multivalued Hamiltonian 
and asymptotically fiat black hole solutions that are stable against Hawking evaporation. 

We take the only bulk contributions to the action to be the Einstein-Hilbert term and 



the four- dimensional Euler density. In D spacetime dimensions, the action then reads [|TT 

+ (boundary terms) , (1.1) 
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where k is the dimensional gravitational constant and A is the single Lovelock parameter 
For D > 5, the four-dimensional Euler density contributes to the equations of motion, and 
we obtain a one-parameter family of generalizations of Einstein's theory. In these theories, 
asymptotically fiat black hole solutions that are stable against Hawking evaporation occur 



only when D = 5 and A > [|TT|], and we shall therefore concentrate on this case. For the 
interest of comparison, we shall also include the limiting case of five-dimensional Einstein 
theory, D = 5 and A = 0. 

We shall formulate the spherically symmetric Hamiltonian theory with thermodynami- 
cally motivated boundary conditions similar to those introduced in Ref . [ET| . On a classical 



solution, one end of the spacelike hypersurfaces will be at the bifurcation three-sphere of a 
nondegenerate Killing horizon, and the other end will be on a timelike hypersurface in an 
exterior region of the spacetime. We shall refer to the two ends respectively as the "left" end 
and the "right" end, as motivated by the Penrose diagram in which our classical solutions 
are embedded in the right-hand-side asymptotically fiat region [|TT|: on a solution, one can 



think of the left end as the inner one and the right end as the outer one. At the left end we 
fix the rate at which the spacelike hypersurfaces are boosted with respect to the coordinate 
time, and at the right end we fix the intrinsic metric on the timelike hypersurface. 

For A > 0, the super-Hamiltonian H turns out to be a multivalued function of the canoni- 
cal variables. However, our boundary conditions are sufficient to uniquely determine H near 
the left end of the spacelike hypersurfaces, and this solution for H can then be uniquely 
extended to the full spacelike hypersurfaces by continuity. Our boundary conditions at the 
horizon thus eliminate the difficulties due to the multivaluedness of the super-Hamiltonian. 

We shall find that the theory admits a natural generalization of the canonical trans- 



formation of Refs. ||20|J21[] . The constraints become exceedingly simple, and a Hamiltonian 
reduction leads again to a single canonical pair of unconstrained degrees of freedom. On 
a classical solution, one member of the pair is the mass parameter, and its conjugate mo- 
mentum is the difference of the Killing times at the left and right ends of the spacelike 
hypersurfaces. 

After taking the curvature radius at the right end of the hypersurfaces to be time- 



independent, we quantize the reduced theory by Hamiltonian methods. Following Ref. ^1 
we analytically continue the time evolution operator to imaginary time and take the trace, 
and interpret the resulting object as the partition function of a thermodynamical canonical 
ensemble. This ensemble describes black hole spacetimes in a spherical "box" whose size 
and boundary temperature are fixed. 

In the special case of Einstein's theory, A = 0, we find that the thermodynamical proper- 
ties of the system are highly similar to those of the corresponding system in four dimensions 



1 2T| , |36|j37| . For high boundary temperatures, the partition function is dominated by a black 
hole that fills most of the box. For low boundary temperatures, on the other hand, there is 
no dominant classical solution. One can argue that the behavior of the partition function 



suggests a topological phase transition from a black hole to "hot flat space" [p6ip7| . 



^We have set c = h = 1: the gravitational constant k has the dimension of (length)^ ^, and the 
Lovelock parameter A has the dimension of (length)'^. 
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For A > 0, the partition function displays several qualitatively different regions depending 
on the relative magnitudes of the box, the temperature, and A. In the high temperature 
limit, with the other two parameters fixed, the partition function is again dominated by a 
classical black hole solution that fills most of the box. In the low temperature limit, with 
the other two parameters fixed, the partition function is now also dominated by a black hole 
solution: this black hole is small compared with the box, and it has no analogue in Einstein's 
theory. However, if A is small compared with the size of the box, the existence of the new 
dominating solution at low temperatures only has a minor effect on the behavior of the 
partition function. In this sense, we can say that the qualitative thermodynamical behavior 
of the pure Einstein system is stable against the addition of the Lovelock parameter. 

When the size of the box is taken to infinity, the partition function does not have a 
well-defined limit, neither for A = nor for A > 0. For A = this is not surprising: 
just as in four dimensions, it reflects the fact that a Schwarzschild hole in asymptotically 
flat space is not stable against Hawking evaporation For A > 0, on the other 

hand, the theory does admit asymptotically fiat black hole solutions that are stable against 
Hawking evaporation |[Tl|, and one might therefore have expected the infinite box limit to 
exist. The reason why this is not the case becomes apparent when one tries to repeat our 
analysis with boundary conditions that replace the right-hand-side timelike boundary by 
an asymptotically flat infinity. The classical reduction and the construction of a quantum 
theory proceed without difficulty, but the effective Euclidean action of the system turns 
out to be unbounded below, and the formal integral expression for the partition function is 
divergent. The effective action has a local minimum, corresponding to the black hole that is 
stable against Hawking evaporation |^8|], but this is not sufficient to ensure the existence of 
the full canonical ensemble. Another system with locally thermodynamically stable classical 
solutions but no well-defined canonical ensemble is four- dimensional Einstein-Maxwell theory 
with fixed charge in asymptotically flat space ||2^ . 

The rest of the paper is as follows. In section |I| we introduce the notation and present 
the Hamiltonian formulation in the metric variables. In section |T| we perform a canonical 
transformation to Kuchaf-type variables and reduce the theory to its unconstrained Hamil- 
tonian degrees of freedom. The reduced theory is quantized and the partition function 
constructed in section and the thermodynamics is analyzed in section 0. Section jVl 
presents a brief summary and discussion. Appendix ^ reviews briefly the classical solutions 
to the theory. The classical and quantum mechanical analyses under asymptotically flat 
boundary conditions are outlined in appendix ^. 

For the rest of paper, we shall use Planckian units that have been tailored for numer- 
ical convenience: we set c = h = 1 and k, = Gvr^. All dimensionful constants have thus 
disappeared. 



II. METRIC FORMULATION 

In this section we introduce the model and present the Hamiltonian formulation in the 
metric variables. 

We begin with the general five- dimensional spherically symmetric Arnowitt-Deser-Misner 
(ADM) metric. 
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-N^df + K^{dr + N'dtf + R^dVt 



(2.1) 



Here dVL\ is the metric on the unit three-sphere, and A^, A^'', A, and R depend on the 
coordinates t and r only. The coordinate r has the range < r < 1; this is convenient in 
view of our boundary conditions, which will make the radial proper distance on the constant 
t hypersurfaces finite. Unless otherwise stated, we assume both the spatial metric and the 
spacetime metric to be nondegenerate. In particular, we take A, i?, and A^ to be positive. 

Inserting the metric ( p.l|) in the Lovelock action ( |1.1|) with Z) = 5, integrating over the 
three-sphere, and dropping a total derivative, we recover the action 



where 



A - (A^^A)'l {R - N-'R'^ , _2 



dt j dr C 



(2.2) 



{R - N'-R' 
N 



+NAR 



N 

1 r 



kR-\ 





/R'Y iR-WR'f 
A + 3iV^ 




The overdot and the prime denote respectively 8/ {dt) and d/{dr) 
conforming to the notation of Ref. 
of motion obtained from local variations of 5"^ 



(2.3) 



We have written A = -^A, 



We have verified that the Lagrangian equations 
]2|) are equivalent to the full spherically 
symmetric Lovelock equations derived from the action ( |1.1| ). The reduction of the 

action by spherical symmetry is therefore consistent with the equations of motion, and we 
can take Sj^ ( |2.2| ) as the starting point of the dynamical analysis. We shall address the 
boundary conditions and boundary terms within the Hamiltonian formulation below. 

For the reasons discussed in section |I|, we take A > 0. For presentational simplicity, we 
shall assume A > until explicitly stated otherwise in section ^ In the limiting case of 
five-dimensional Einstein gravity, A = 0, the analysis would proceed in an entirely analogous 
manner, with the obvious technical simplifications. 
The Hamiltonian form of the action (E.2|) is 



(2.4) 



Ss = J dt dr (PaA + PrR -NH - N''H, 



where the super-Hamiltonian constraint and the supermomentum constraint are given re- 
spectively by 



H = y{PR + y 



Hr 




(2.5a) 
(2.5b) 
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The quantity y is determined in terms of the canonical variables by the cubic equation 



= \\y^ + y}R^ + \ 




+ Pa 



(2.6) 



Note that the form of the supermomentum (|2.5b|) is completely determined by the fact that 
it must generate spatial diffeomorphisms in all the canonical variables, together with the 
observation that R and Pa are spatial scalars whereas A and Pr are spatial densities . 

Depending on the values of A, R, and Pa, the cubic ( p. 61) can have up to three real 
solutions for y. The super-Hamiltonian is therefore a potentially multivalued function of 
the canonical variables. Such multivaluedness arises generically in Lovelock theories, owing 
to the presence of kinetic terms higher than quadratic in the velocities in the Lagrangian 
density in our case, the highest kinetic terms in ( p.3|) are quartic in the velocities. 

We shall address this phenomenon in more detail below. The geometrical meaning of y is 
revealed by observing that the equation of motion obtained by varying S'e ( |2.4| ) with respect 
to Pr reads 



R = Ny + N'R! 



(2.7) 



On a classical solution, y is therefore uniquely determined by the embedding of the spacelike 
hypersurface in the spacetime. Conversely, when multiple real solutions to ( |2.6| ) exist, it can 
be verified that they generically lead to different spacetimes. 

Let us turn to the boundary conditions. From the Lovelock generalization of Birkhoff 's 
theorem it follows that the local properties of the classical solutions are completely 
characterized by a discrete binary parameter and a continuous, mass-like parameter. The 
general solution is shown in curvature coordinates in appendix 0. We wish to concentrate 
on the black hole solutions, whose global structure is similar to that of Kruskal manifold 
Tl| . We further wish to attach the left end of our spacelike hypersurfaces to the bifurcation 



three-sphere, and to prescribe there the rate at which the hypersurfaces are boosted with 
respect to our coordinate time. The right end of the hypersurfaces will then be in the right- 
hand-side exterior region, and we wish to prescribe the metric on the timelike hypersurface 
that this end traces. We must now specify boundary conditions and boundary terms that 
achieve this. 

Consider first the left end of the hypersurfaces. Following the analogous treatment in 
Refs. 



2T|J2|J2|J2§, we adopt at r ^ the falloff 



(2.8a) 
(2.8b) 
(2.8c) 
(2.8d) 
(2.8e) 
(2.8f) 

where Aq and Rq are positive, and Ni > 0. 0(r") stands for a term that is bounded at r ^ 
by a constant times r", and whose derivatives fall off accordingly. As in Refs. pT|,pB|,|2^,pB| , 



A(t,r) 


= Ao(t) + 0(r2) , 


R{t,r) 


= Po(t) + i?2(t)r' + 0(r^) 




= 0{r') , 


Pnit^r) 


= 0(r) , 


N{t,r) 


= Ni{t)r + 0{r^) , 


N''{t,r) 


= Nl{t)r + Oir^) , 
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these conditions guarantee that the classical solutions have a bifurcate horizon, they put the 
left end of the spacelike hypersurfaces at the bifurcation three-sphere, and they are consistent 
with the constraints and preserved by the Hamiltonian evolution. They also ensure that the 
cubic ( p. 61) has a unique real solution for y near r = 0. On a classical solution, the future 
unit normal vector n'^{t) to the spacelike hypersurfaces at r = then evolves according to 



cosh 



t2 



AQ\t)Ni{t)dt 



(2.9) 



Next, consider the boundary conditions in the variational principle. At r = 0, we follow 
Refs. [^^PHp6| and make Ni/Aq a prescribed function of t. By (|2.9| ), this means fixing 
the rate at which the constant t hypersurfaces are boosted at r = 0. At r = 1, we make R 
and —gtt = N"^ — (AN^)'^ prescribed positive- valued functions of t. This means fixing the 
intrinsic metric on the three-surface r = 1, and in particular fixing this metric to be timelike. 



Finally, we need the boundary terms to be added to the bulk action (|2.4| ). As in Ref. 
I2TII, it can be verified that the appropriate term at r = is 



dtRo {IRI + X) (A^i/Ao 



and the appropriate term at r = 1 is the integral over t of 
NA-^R^R' - N' APa - \R{R^ + A) In 



(2.10) 



N -AN"^ 





R{R - N'-R') 



Ar2 



AA^^A 
3A^ 



R' 



m - {AN'-y 



+ 



N^'iR') 



(2.11) 



We have therefore arrived at a variational principle with the desired boundary conditions. 
The Lovelock generalization of Birkhoff 's theorem guarantees that classical solutions exist, 
and makes possible a complete description of the solutions. One should note, however, that 
our Hamiltonian action does not directly reflect the split of the variables at r = 1 into the 
dynamical degrees of freedom versus the boundary data. With the data on the timelike 
boundary, the evolution no longer forms a hyperbolic system. 



III. CANONICAL TRANSFORMATION AND HAMILTONIAN REDUCTION 

In this section we simplify the constraints by a canonical transformation and reduce 
the theory to unconstrained Hamiltonian variables. The treatment will closely follow Refs. 

To begin, suppose that we are given the canonical data (A, R, P\, Pr) on a spacelike 
hypersurface embedded in a classical solution. We wish to reconstruct from this data both 
the spacetime and the location of the hypersurface in the spacetime. 

As we have noted above, the embedding of the hypersurface in the classical solution 
defines a unambiguous value of y: from the equation of motion (p^.T]), one finds that this 
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value is ytme := N~^{R — N^R'). To reconstruct ytme from the canonical data, one needs to 
solve the cubic ( p.6|) , which may have up to three real solutions. Near r = 0, the falloff ( p.8|) 
guarantees that the cubic has a unique real solution, and this solution must therefore be equal 
to i/true- As r increases, two spurious real solutions may appear, but it is straightforward to 
verify that neither of the spurious real solutions can ever be equal to ytme- Therefore, ytrue 
is recovered from (|2.6| ) by choosing the unique real root near r = and following this root 
by continuity to all r. We note that, generically, neither of the spurious roots for y satisfies 
the constraint H = 0. 

After y = ytme has been recovered, the reconstruction proceeds in full analogy with that 
in Ref. ||2^. The function F appearing in the metric ([Al|) is given by 




y' , (3.1) 



and from ([A2D one finds for the mass the expression 

M=^R\l-F) + lX{l-Ff . (3.2) 



Finally, one finds 



T' = ^ , (3.3) 



which specifies the location of the hypersurface up to translations in the Killing time. This 
completes the reconstruction. 

Next, we wish to promote the reconstruction equations into a canonical transformation, 
valid even when the equations of motion do not hold. Provided we stay within a sufficiently 
narrow neighborhood of the classical solutions, y is again uniquely recovered as a function 
of the canonical data by taking the unique real root of ( |2.6| ) near r = and continuously 
following this root as r increases. Computing the Poisson bracket between M and T' suggests 
that —T' could serve as the momentum conjugate to M; if this holds, the new momentum 
conjugate to R := i? is fixed by the fact that the supermomentum constraint generates spatial 
diffeomorphisms in all the variables and must thus read PmM' + PrR' . These considerations 
suggest the transformation 

M := lR\l - F) + fX{l - Ff , (3.4a) 

Pm--=-^ , (3.4b) 
R :=R , (3.4c) 
PR:=F--'[A~^R'Hr-yH) , (3.4d) 

with F given by ( |3.1| ). We now need to examine whether this transformation is indeed 
canonical. 

To proceed, we arrange the difference of the integrands in the Liouville forms as 
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Pa5A+Pr5R - PmSM - PrSR 
= 6 



+ 



APa - XyA-\R'f + \R' (r^ + A) In 

R' + yA 



XyA-^R! -UR' + 



R'-yA 




(3.5) 



Both terms on the right-hand side of (3^) are well defined. Upon integration from r = 
to r = 1, the second term only produces contributions from the two ends. The contribution 
from r = vanishes because of the falloff (|2.8|) . The contribution from r = 1 vanishes if 6R 
vanishes there. As 6 should in the context of the Liouville form be understood as a time 
derivative, this happens when the boundary conditions fix R to be independent of t at r = 1. 
If this is the case, we see that the difference of the Liouville forms is an exact form. 



dr{PA6A + Pr6R) - dr {Pm6M + Pr5R) 
Jo Jo 



APa - XyA-\Rf + \R! (r^ + A) In 




R' + yA 



R' - yA 



(3.6) 



and the transformation is canonical. 

If, on the other hand, the boundary conditions fix R to be explicitly t-dependent at 
r = 1, one cannot similarly argue that 5R would vanish at r = l.Q As mentioned at the 
end of section ||, the canonical variables at r = 1 do not cleanly split into "independent" 
degrees of freedom versus boundary data, and it us unclear to us what the proper attitude 
here should be. We shall, nevertheless, proceed to regard the transformation as canonical 
even when R is explicitly t-dependent at r = 1: as in Refs. ||2ll , |23| , |26| , it will be seen that 
no apparent inconsistency will result. From the viewpoint of thermodynamics, the case of 
principal interest will in any case be the one where R is independent of t at r = 1. 

By construction, our transformation is well defined in a sufficiently narrow neighborhood 
of the classical solutions. It also has a unique inverse. Equations ( |3.4a| ) and (|3.4c|) , together 
with the falloff implied by ( ^.8[ ), determine F uniquely in terms of M and R. Equations 
( p.l| ) and (|3.4b| ), together with the fact that A is by assumption positive, then determine A 
and y. P\ is obtained from (p.6| ), and Pr finally from ( |3.4d|) . 

To obtain the action in the new variables, we note that the constraint terms can be 
written as 



NH + N'^Hr = N^M' + iV'PR 



where 



N 



M 



-NP-^A'^R! - N''F-^Ay 



= Ny + N^R! 



(3.7) 



(3.8a) 
(3.8b) 



^This appears to have gone unmentioned in Refs. |21,p3,2 
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This suggests that one could take A^-^ and A^*^ as the new independent Lagrange multiphers 
in the action. Examining the falloff at r = reveals, however, that fixing A^*'^ at r = to 
a value that is independent of the canonical variables is not equivalent to fixing A'^iAg ^ to 
a value that is independent of the canonical variables. This difficulty can be remedied by 
redefining the Lagrange multipliers near r = as in Ref. |2^, and the appropriate boundary 
terms at r = and r = 1 can then be constructed as in Refs. [^^^ . After these steps, 



the constraints can be eliminated by a Hamiltonian reduction as in Refs. pmpT] , pB[ , and one 
recovers a reduced theory in a true Hamiltonian form. The steps follow the cited references 
so closely that we shall here omit the detail and proceed directly to the reduced action. 
The reduced action reads 



•S'red = J (P"^ ~ h) 



(3.9) 



The coordinate m arises from the unreduced theory as the r-independent value that M takes 
when the constraint M' = holds. The momentum p is related to the unreduced variables 
by 



Jo 

The Hamiltonian h is given by 

h = -NoRh {iRl + A 



drP, 



M 



(3.10) 



Q2F + 52 + i51n 



'Q^F + B^ - B 



lQ'^? + B^ + B, 



3/2 



where 



No:=Ari/Ao , 



(3.11) 



(3.12) 



Rh 



2m- iA 



(3.13) 




(3.14) 



and B and are respectively the values of R and —gu at r = 1. B, Q^, and Nq are 
considered prescribed functions of t, satisfying i? > 0, > 0, and Nq > 0. The range of m 
is |A < m < ii?^ + iA, corresponding to Q < Rh < B, and the range of p is the full real 
axis. Rh is equal to Rq, and on a classical solution it equals the horizon radius. 

The equation of motion for m implies that m is independent of t: the value of m is 
simply the mass parameter of the classical solution. The equation of motion for p refiects 
the fact that, by (|3.3| ), ( p.4b|) , and ( |3.1CI| ), p is equal to the difference in the Killing times at 
the two ends of the spacelike hypersurface. 
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IV. QUANTIZATION AND THE PARTITION FUNCTION 



In this section we quantize the reduced Hamiltonian theory and obtain a partition func- 
tion as the trace of the analytically continued time evolution operator. 

From now on, we take the boundary radius independent of time, B = 0. We also subtract 
from the Hamiltonian the value that the terms arising from r = 1 would take on flat 

spacetime. This subtraction does not affect the equations of motion, but it does renormalize 
the value of the action: it is analogous to subtracting the Kq term in Einstein's theory 
39| , |40[| . Writing Q := \fQ^ > 0, the new Hamiltonian is given by 



h = g(l - v^) [b^ + |A (l - v^) (2 + v^)] - UoRh {iRl + A 



(4.1) 



The first of the two terms in ([4.1D is the Lovelock analogue of the quasilocal energy of Brown 



and York pl| -^. The second term arises from the bifurcation three-sphere, and it will give 
rise to the black hole entropy. 

We take the wave functions to 



Quantization proceeds exactly as in Refs. ||2T|j2^j26 



be functions of the configuration variable m, with < m < ^B"^ + jX, and we introduce 
an inner product with some smooth and slowly varying weight factor. The Hamiltonian 
operator is taken to act by multiplication by the function h ( [4.1| ), iplm) h(m)'?/'(m), and 
the unitary time evolution operator is easily found. We then analytically continue the argu- 
ments of the time evolution operator to imaginary values: we set / Qdt = —ijS, interpreting 
/3 > as the inverse temperature at the boundary, and / NqcH = —2711, motivated by the 
regularity of the classical Euclidean solutions. The trace of the analytically continued time 
evolution operator is divergent, but we can argue as in Refs. pT| , |2^p6| that an acceptable 
renormalization is achieved by introducing a suitable regularization, dividing by the trace 
of the regularized identity operator, and finally eliminating the regulator. In this fashion, 
we obtain for the renormalized trace the manifestly well-defined expression 



Z{P-B-X) 



fidx 



where the effective action is given by^ 



= /3B^ (l - Vf' 



A 



with 



fidx exp{—I^:] 



1- v^U2 + v^' 



2nB'x Ux' + A 




(4.2) 



(4.3) 



(4.4) 



We have introduced the integration variable x = Rh/B, and the smooth and slowly varying 
positive function fi{x) arose from the choice of the inner product. 



^This effective action has been obtained previously [E3 by the Euclidean Hamiltonian reduction 



method of Ref. |37|. 
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We now interpret the object Z{(3; B; A) (|4.2| ) as the partition function of a thermodynam- 
ical canonical ensemble describing black holes in a spherical box with curvature radius B 
and and inverse boundary temperature /3. The thermodynamical properties of this ensemble 
will be analyzed in the next section. 



V. THERMODYNAMICS IN THE CANONICAL ENSEMBLE 



As noted above, the partition function Z{/3; B; A) ( [4 .21) is manifestly well defined. Fur- 
ther, the form of the integral in (|4.2|) guarantees that the (constant volume) heat capacity, 
C = P'^(d^{\n Z)/dP^), is always positive (see, for example, section IV of Ref. ||26l), and 
that the ensemble has a well-defined density of states p6|- p8| , ^ . These properties support 
the interpretation of the partition function in terms of a genuine thermodynamical equilib- 
rium ensemble, in spite of the fact that we arrived at the partition function via an analytic 
continuation and not via direct statistical mechanics arguments. 

To proceed, we shall estimate the integral in ( |4.2| ) by the saddle point approximation. 
We shall throughout assume fi{x) to be so slowly varying that its precise form will not 
affect the saddle point analysis. We shall also assume that the action is sufficiently rapidly 
varying to make the saddle point approximation is justified, without attempting to explicitly 
state the necessary conditions; typically, it will be throughout assumed that the system is 
"macroscopic," B ^ 1. 

The critical points of are at the roots of the equation 



X 



(5.1) 



The critical points give precisely the Lorentzian black hole solutions whose Hawking tem- 



perature at the boundary, calculated in the usual way from the surface gravity and the 
blueshift factor, is equal to (3. The mass of the hole is m = ^B'^x'^ + |A, and the value of 
at a critical point equals the Euclidean action of the corresponding Euclideanized black 
hole solution. Whenever the partition function is dominated by a critical point, we recover 
for the thermal energy expectation value and the entropy the results 



S 



B^ (l - Vf' 



A 



1 - v^) (2 + v^' 



d{\nZ) 

1 - /^^) (In^) ^ 27ri?3x Ux' + ^ 1 = 2nR, (ii?^ + A 



(5.2a) 
(5.2b) 



where x and F are evaluated at the critical point. The expression ( p.2b| ) for the entropy 



agrees with the result first obtained by Euclidean methods [|Tl 



We can now extract physical information by analyzing the critical point structure of 
in various limits of interest in the three parameters A, B, and (3. 

As a preliminary, consider the case A = 0, in which our Lovelock theory reduces to 
Einstein's theory. Although we have for presentational simplicity assumed A > 0, it is easy 
to see that the partition function for Einstein's theory is correctly recovered by taking the 
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limit A — > in equations ( [4.2|) -( ^!4D . In particular, (|4.4| ) reduces to F = 1 — x^, and the 
critical point equation ( |5.1| ) reduces to 



27r5 



X 



Vl -x2 . (5.3) 



X 



The condition for critical points to exist is f3 < tiB, and the critical points are then at 
x± := 2~^/^(^l ± VI ~ Ti^'^ B^'^ 13'^^ . When the critical points are distinct, x+ is a 
local minimum and x_ a local maximum. When < ^tt'^B^, the partition function gets 
its dominant contribution from the global minimum at x = x+. When > ^tt^B"^, on 
the other hand, the partition function gets its dominant contribution from the vicinity of 
the global minimum at x = 0. The limiting case = ^ti'^B^ represents a phase transition 
where the dominant contribution shifts from x = x+ to x = as /3 increases. When the 
saddle point dominates, the thermal energy and entropy (p. 21) take the form 



{E) ^ B' i^l- ,Jl-xlj , (5.4a) 
S ^ IttB^x'I = IttRI , (5.4b) 
and the relation between the thermal energy and the mass can be written as 

m^{E)-^^ . (5.5) 

Equation (|5.5|) displays explicitly how the mass gets a contribution both from the thermal en- 
ergy and from the gravitational binding energy associated with the thermal energy. Expect- 
edly, the situation is closely similar to that in four- dimensional Einstein theory [p6| , p7| , ^ . 



We now turn to the case A > 0, in which J* always has at least one critical point. 

Consider first the limit of small A with fixed B and f3. The situation differs from that 
in the case A = only in that there is now one new critical point, a local minimum, at 
X = 2'kXB-^(3-^ + 0{X'^). At the new critical point, h = IX/S + OiX"^). Therefore, as A 0, 
the partition function smoothly approaches that of Einstein's theory. In particular, when 
P"^ < jTi'^B'^, it would be straightforward to compute the first order correction in A to the 



thermal energy and the entropy ( |5.4| ), assuming that the corrections to the saddle point 
approximation are small. 

Consider next the small (3 limit with fixed B and A. There is only one critical point, 
at X = 1 — |7r~^(i?^ + A) + 0(/?^), and this critical point is the global minimum of J*. 
One can think of this critical point as the counterpart of the larger of the two critical points 
of the case A = 0: the black hole fills almost all of the box. The disappearance of the 
smaller critical point of the case A = is related to the fact that, for fixed A, the Hawking 
temperature of the Lovelock hole in asymptotically fiat space is bounded below by |7r~^A~^/^ 
pj| . If the saddle point approximation to the partition function remains good, the thermal 



energy and the entropy are given by 

^Bf3 + 0{(3^) , (5.6a) 




8n^B^ 



0{p^) . (5.6b) 
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Consider next the large j3 limit with fixed B and A. There is again only one critical point, 



at X = 27rAfi-i/5-iFy^ + 0(/3-3)^ where Fq := 1 + B'^X-^ [1-^1 + X^B-^j. This critical 

point is the global minimum, and it has no counterpart in Einstein's theory: it corresponds 
to a small, "purely Lovelock," black hole. If the saddle point approximation to the partition 
function remains good, the thermal energy and the entropy are easily read off from (|5.2|) as 



{E)^B^i-?r 



+ 0(/3-') , (5.7a) 



_J^ + 0(r') • (5.7b) 

Consider then the large A limit with fixed B and (3. There is again only one critical 
point, at a; = 1 — |7r^^/5^A^^ + 0(A^^), and this critical point is the global minimum. The 
hole is again "purely Lovelock," but it now fills almost all of the box. 

Finally, consider the large B limit with fixed A and (3. One critical point is at x = 
1 — |7r~^i?~^/3^+0(i?~^). This critical point is the global minimum, and it can be regarded as 
the counterpart of the larger of the two critical points of the case A = 0. \i (3 > 47rA^/^, there 



are in addition two other critical points, at x = ^tt ^B ^1 ± y 1 — IGvr^A/? +0{B ^). 

The fact that the two small critical points exist only for (3 > AnX^^"^ is related to the above- 
mentioned phenomenon that the Hawking temperature of our Lovelock hole in asymptoti- 
cally fiat space is bounded below by ^n^^X^^^"^ ||Tl|. If the saddle point approximation is 



good, the thermal energy and the entropy are obtained by replacing both 0-terms in ( ^.6|) 
byO(S-i). 

We therefore see that for A > 0, the partition function is always dominated by a black 
hole solution in the limits that we have considered. In the high temperature limit and in 
the large box limit, the situation is very similar to that for A = 0, in that the dominating 
black hole solution fills most of the box. In the low temperature limit, on the other hand, 
the Lovelock theory does exhibit a dominating black hole solution where none existed in 
the case A = 0. For a macroscopic box and A <C -B^, however, the presence of the new 
dominating solution does not appear to make the thermodynamical behavior qualitatively 
different from that in the case A = 0. One can read these results as evidence for stability 
of the qualitative thermodynamical behavior of Einstein's theory upon the addition of the 
Lovelock parameter. 

For A > 0, the critical point structure of /* is entirely determined by two parameters, 
which can be conveniently taken to be XB~^ and (3B~^. Numerical experimentation suggests 
that there are never more than three critical points. When XB~^ is fixed and sufficiently 
large, there is only one critical point for any (3B~^: this critical point is the global minimum, 
and it migrates smoothly from the large Lovelock hole to the small Lovelock hole as f3B~^ 
increases. When XB~^ is fixed and sufficiently small, on the other hand, the transition 
between the unique minima for small and large j3B~^ takes place via a phase in which /* 
has three critical points, a maximum surrounded by two minima: as (3B-^ tends to zero 
(infinity), only the larger (smaller) of the two minima prevails. When three critical points 
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exist, one can find regions of tfie parameter space wfiere tlie global minimum is at either 
of the two local minima. We have, however, not attempted to corroborate these numerical 
experiments analytically. 

It should be emphasized that the partition function has no well-defined limit as B ^ oo 
with fixed (3 and A, neither for A = nor for A > 0. As with Einstein's theory in four 
dimensions ^J3^, this reflects the fact that the thermodynamical canonical ensemble is 
not well defined in asymptotically fiat space. We shall give a more detailed comparison of 
the boxed Lovelock theory to Lovelock theory in asymptotically fiat space in appendix R 



VI. SUMMARY AND DISCUSSION 

In this paper we have investigated the Hamiltonian dynamics and thermodynamics of 
five-dimensional spherically symmetric Lovelock theories in which the only contributions 
to the Lagrangian density are the Einstein-Hilbert term and the four-dimensional Euler 
density. We adopted boundary conditions that enforce every classical solution to be part 
of the exterior region of a black hole, with the spacelike hypersurfaces extending from the 
horizon bifurcation three-sphere to a timelike boundary with fixed intrinsic metric. We 
simplified the constraints by a canonical transformation that generalizes the one introduced 
by Kuchaf in four- dimensional spherically symmetric Einstein theory, and we reduced the 
theory classically to its true dynamical degrees of freedom. 

After Hamiltonian quantization, we interpreted the trace of the analytically continued 
time evolution operator as the partition function of a thermodynamical canonical ensemble, 
describing black holes in a spherical box whose size and boundary temperature are fixed. In 
the special case where the Lovelock parameter A vanishes and the theory reduces to Einstein's 
theory, we found that the thermodynamics is highly similar to that of the corresponding 
system in four- dimensional Einstein theory: in particular, for high boundary temperatures 
the partition function is dominated by a classical black hole solution that fills most of the 
box. When A > 0, the situation was more versatile. In the high temperature limit, with A 
and the box size fixed, the partition function is again dominated by a black hole that fills 
most of the box. In the low temperature limit, on the other hand, the partition function 
is now also dominated by a black hole solution; this black hole is small, and it has has no 
analogue in Einstein theory. Nevertheless, if A is small compared with the size of the box, 
the new dominating solution has little qualitative effect on the thermodynamical properties. 
In this sense, the qualitative thermodynamical behavior of the Einstein system is stable 
upon the addition of the Lovelock parameter. 

When the box size is taken to infinity, we found that the partition function has no well- 
defined limit, neither for A = nor for A > 0. While this is not surprising for Einstein's 
theory, in view of the similar phenomenon in four dimensions [^,^], one might have hoped 
the theory with A > to fare better on the grounds that this theory admits asymptotically 



fiat black hole solutions that are stable against Hawking evaporation [|ri[|. However, even 



though a classical solution that dominates a well-defined partition function must be stable 



against Hawking evaporation our Lovelock theory in asymptotically flat space provides 
an example where the mere existence of such a locally stable classical solution does not 
imply the existence of well-defined canonical ensemble. Another such example occurs in 
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four- dimensional Einstein-Maxwell theory in asymptotically fiat space [p4 



In the classical theory with A > 0, we saw that the super-Hamiltonian emerges as a 
multivalued function of the canonical variables, as is generically the case in Lovelock theories 
]3J|]. Nevertheless, our thermodynamically motivated boundary conditions were sufficient 
to uniquely specify the super-Hamiltonian near the horizon, and the uniqueness could be 
extended to the full spacelike hypersurfaces by continuity. Another boundary condition that 
would uniquely specify the super-Hamiltonian in this fashion is the asymptotically flat falloff 
( pi|) discussed in appendix |B[ However, one expects there to exist boundary conditions of 
interest for which such uniqueness does not occur, and in such cases one would need to seek 
other criteria for specifying the super-Hamiltonian. If one regards the Lovelock theory as 
a perturbation to Einstein's theory, or as a toy model for semiclassical gravity with back- 
reaction, one possible criterion of this kind is perturbative expandability of the solutions 
in A [H. 

In conclusion, our results provide evidence for robustness of the classical Hamiltonian 
structure and the qualitative thermodynamical structure of spherically symmetric Einstein 
gravity in five dimensions upon the addition of the four- dimensional Euler density in the 
action. To put this conclusion in proper perspective, one should remember that both our 
particular Lovelock theory and our boundary conditions were hand-picked so that the global 
aspects of the problem remained virtually identical to those in pure Einstein gravity. It is 
tempting to think that this may exemplify a more general connection between the global 
properties of the (space of) classical solutions and the qualitative behavior of thermody- 
namical equilibrium ensembles: one might conjecture that whenever the global properties 
of a Lovelock theory are sufficiently similar to those of Einstein's theory, then also the equi- 
librium thermodynamics, with finite or infinite boundary conditions, will be qualitatively 
similar to that in Einstein's theory. Another example supporting such a conjecture is pro- 
vided by the asymptotically anti-de Sitter Lovelock theories of Ref. [^, which include as 
a special case Einstein gravity in three and four dimensions with a negative cosmological 
constant. However, to give the conjecture a more substantial meaning, one would need a 
more systematic understanding of the possible global structures that the various Lovelock 
theories may have. 
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APPENDIX A: LOVELOCK BLACK HOLE 



In this appendix we briefly review the classical solutions to our Lovelock theory. 
By the Lovelock generalization of Birkhoff's theorem ||^, the general solution to the 
theory for A 7^ can be written in the local curvature coordinates (T, R) as 



ds^ = -FdT^ + F-^dR^ + R^dnl , (Al) 

where 

(A2) 




with A = |A. The coordinates T and R are respectively the Killing time, whose constant 
value hypersurfaces are spacelike or timelike depending on the sign of F, and the curvature 
radius of the three-sphere. For A = 0, the general solution is the five- dimensional Schwarz- 
schild solution, obtained from (|A2|) with the lower sign as the limit A 0. In addition to A, 



the only local parameter in the solution is the mass-like quantity cu. 

For the reasons discussed in section |, we take A > 0. The solution then describes a 
black hole in asymptotically fiat space provided we take > |A and, for A > 0, choose 
the lower sign in (|A2|) [|lT|. The curvature coordinates are good individually in each region 
not containing horizons: the horizons are nondegenerate, and the Penrose diagram of the 
conventional maximal analytic extension is similar to that of Kruskal manifold. The horizon 
radius is Rh = — |A. In our units, the ADM mass is M = 



APPENDIX B: ASYMPTOTICALLY FLAT INFINITY 



In this appendix we adapt the analysis of the main text to boundary conditions that re- 
place the timelike boundary by an asymptotically flat infinity. We shall see that quantization 
along the lines of section |V| will not lead to a well-defined canonical ensemble. 

In the metric theory of section ^ we let r take the range < r < 00. At r ^ 00, we 
introduce the falloff 



A(t,r 
R{t,r 
PAit,r 

N{t,r 
N'{t,r 



1 + M+{t)r-^ 
O 



^oo 1^-1- 



(Bla) 
(Bib) 
(Blc) 
(Bid) 
(Ble) 
(Blf) 



where N+{t) > 0, and e is a parameter that can be chosen arbitrarily in the range < e < 2. 

00 by a constant times r'^, and whose 



fr"*) denotes a term that is bounded at r 
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derivatives fall off accordingly. It is straightforward to verify that this falloff makes the 
coordinates asymptotic to hyperspherical coordinates in Minkowski space, it is consistent 
with the constraints, and it is preserved by the time evolution. iV+(t) gives the rate at which 
the asymptotic Minkowski time advances with respect to the coordinate time t. When the 
equations of motion hold, M_|_(t) is independent of t, and its value is the ADM mass. 
The total action takes the form S = Sj] + Sqj], where 

5as = JdtRo {iRl + a) (NjAo) - J dt , (B2) 



and Sy, is as in ( |2.4| ) except that the upper limit of the r-integral is replaced by infinity. The 
canonical transformation and Hamiltonian reduction proceed as in section |nT[ The action 
of the reduced theory is as in (|3.9|) , but with the Hamiltonian now given by 



\i = N+m-HoRh{\Rl + \) . (B3) 

The configuration variable m has the range m > |A, while the range of p is the full real 
axis. Although in the main text we assumed for presentational simplicity A > 0, it is easily 
verified that the above derivation extends to the case A = 0, and the reduced Hamiltonian 
O) is valid for all A > 0. 

Quantization can now proceed as in the main text. We analytically continue the time 
evolution operator by / H^dt = — 27rz and / N^dt = —i(3, interpreting [3 as the inverse 
temperature at the infinity. For the renormalized trace of the analytically continued time 
evolution operator, we obtain formally 

POO 

Z{f3-\)=M lidRheM-IT) , (B4) 
Jo 

where the effective action is given by 

I^ = \^(Rl + \\)-2'nRn{\Rl + \) . (B5) 

The smooth and slowly varying positive function jl{Rh) arose from the choice of the inner 
product, and AT is a normalization constant, possibly dependent on A but presumably not 
on 13. However, the integral in ( p^) is divergent because tends to — oo at large Rh- Thus, 
the canonical ensemble does not exist under the asymptotically flat boundary conditions, 
neither for A = nor A > 0. In this sense, the asymptotically fiat Lovelock theory is 
thermodynamically no better behaved than asymptotically fiat Einstein theory. 

The critical points of give the (Lorentzian) classical solutions that have the inverse 
Hawking temperature (3 at infinity. For A = there exists exactly one critical point, which is 
a local maximum: this is similar to what happens with Einstein's theory in four dimensions 
36| . For A > 0, the situation is more versatile. Critical points exist when (3'^ > IGtt^A, 



and when the inequality is genuine, there are two critical points. The critical point with the 
smaller (larger) value of Rh is a local minimum (maximum, respectively). The local minimum 
gives the classical solution that was found to be stable against Hawking evaporation in Ref. 



11 1 . While the stability of this solution against Hawking evaporation reflects its being a 



local minimum of PS[, the divergence of the integral in (|B4D demonstrates that this local 
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stability is not sufficient to guarantee the existence of the canonical ensemble. The effects 
of the Lovelock parameter on the asymptotically flat thermodynamics are thus qualitatively 
very similar to those of a fixed charge in asymptotically fiat four-dimensional Einstein- 
Maxwell theory 
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